The Astrophysical Journal, 684:1062-1074, 2008 September 10 

Preprint typeset using LM^X style cmulatcapj v. 05/04/06 



TIDAL DEBRIS IN ELLIPTICAL GALAXIES AS TRACERS OF MERGERS WITH DISKS 
R. Feldmann 1 , L. Mayer 2 and C. M. Carollo 1 

ABSTRACT 

We use a set of high-resolution A-body simulations of binary galaxy mergers to show that the mor- 
phologies of the tidal features that are seen around a large fraction of nearby, massive ellipticals in 
the field, cannot be reproduced by equal-mass dissipationless mergers; rather, they are well explained 
by the accretion of disk-dominated galaxies. In particular, the arm- and looplike morphologies of the 
observed tidal debris can only be produced by the kincmatically cold material of the disk components 
of the accreted galaxies. The tidal features that arise from such "cold-accretion" events onto a massive 
elliptical are visible for significantly longer timescales than the features produced by elliptical-elliptical 
mergers (about 1-2 Gyr vs. a few hundred million years). Mass ratios of the order of 1 : 10 between 
the accreting elliptical and the accreted disk are sufficient to match the brightness of the observed 
debris. Furthermore, stellar population synthesis models and simple order-of-magnitude calculations 
indicate that the colors of the tidal features generated in such minor cold-accretion events are rel- 
atively red, in agreement with the observations. The minor cold-accretion events that explain the 
presence, brightness, and structural and color properties of the tidal debris cause only a modest mass 
and luminosity increase in the accreting massive elliptical. These results, coupled with the relative 
statistical frequencies of disk- and bulge-dominated galaxies in the field, suggest that massive ellipti- 
cals assemble most of their mass well before their tidal debris forms through the accretion of relatively 
little, kinematically cold material rather than in very recent, dissipationless major mergers. 
Subject headings: galaxies: interactions — galaxies: elliptical and lenticular, cD — galaxies: spiral 



1. INTRODUCTION 

Bright elliptical galaxies are the most massive stel- 
lar systems in the Universe and host a large frac- 
tion of the stell ar mass that has been produced sinc e 
the Big Bang (jFukugita et alJ 119981 : iRenzinil 120061 ). 
The boxy isophotes and slow internal rotation of the 
most massive ellipticals argue for dissipationless merg- 
ers playi ng a role in the l ast stage of their assembl y 
history (jNaab et alJ 120061 : iKhochfar fc Burkertl I2005D . 
Semianalytic models indicate that major mergers that 
lead to massive ellipticals are typically not disk-disk 
mer gers but rather elliptical-e lliptical or mixed merg- 
ers ([Khochfar fe Burkertl 1200 3). Another scenario, as 
sugge sted by cosmological SPH simulations (|Naab et al.l 
|2007|) . is the early build-up of ellipticals at high redshift 
with only minor growth in mass due to small accretion 
events below z < 1. The red, stellar tidal arms and loops 
in the outskirts of many nearby bright ellipticals are 
widely believed to be fast-dissolving transient beacons 
of recent dissipa tionless mergers bet ween smaller, equal- 
mass ellipticals (|van Dokkum 2005). 3 Such mergers are 
often called "dry" to highlight the absence of gaseous ma- 
terial (as opposed to "wet", or gas-rich, mergers). Nu- 
merical simulations of elliptical-elliptical mergers indi- 
cate that tidal debris disappears w i thin a few hundred 
million years JCombes et al.1 119951: Ivan Dokkum] 120051 : 
iBell et al.ll2006l) . This would imply a substantial increase 



in the massive el liptical galaxy p opulation in the last 
few billion years (jBell et al.l l2004f ). This picture faces, 
however, several observati onal challenges (|Cimatti et al.l 
120061 IScarTata et al.ll2007fl. 
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3 The merger rates derived in lvan Dokkum (2005) are calculated 
from a close galaxy pair analysis and hence should not suffer from 
a misidentification of the origin of the tidal features. 



Other studies ( Kawata et al.l 120061 ) had hinted at the 
gradual disruption of a tiny satellite galaxy by a massive 
elliptical as a source for red tidal features. Yet, these 
disruption events yield tidal debris too faint to match 
the typic al surface brightne sses of the observed tidal 
feature s (Ivan Dokkuml l2005f). E arly numerical experi- 
ments ([Hernquist fc Quinnlll988f ) showed that sharp tidal 
arms and loops could arise from the tidal disruption of 
dynamically cold stellar disks. These experiments how- 
ever adopted a spherical model for the elliptical galaxy 
and lacked the proper resolution to robustly establish 
the morphology and lifetime of the tidal streams. More- 
over, those early studies were carried out before the cur- 
rent cosmological framework, the cold dark matter model 
with a cosmological constant, was established. 

We use a set of 13 high-resolution iV-body simula- 
tions of galaxy-galaxy mergers between a massive ellipti- 
cal galaxy and an elliptical or disk galaxy companion 
to investigate the origin of tidal features observed in 
bright ellipticals. The companion mass ranges between 
10% and 100% of the mass of the main elliptical (Ta- 
ble [l}. The numerical galaxy models are constructed 
to properly matc h the observed stru ctural properties 
of both elliptical (|Bender et al.l 119921 ) and disk galax- 
ies (|Mo et al.lll998l ): in particular, the main massive el- 
liptical has a stellar mass of 2.7 x 10 11 M Q , an effective 
radius of 9 kpc, and a stellar mass-to- light ratio of 4 in 
the i?-band. These parameters ensure that the model 
galaxy li es on the observed fundamental plane of e l- 
lipticals (jDiorgovski fc Davisl[l987l : iDressler et al.lll987h . 
The stellar component of the galaxies is embedded in 
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a cold dark matter halo with properties der ived from 
cosmological simulations ([Navarro et alJ fl.997). The ini- 
tial orbits in the mergers are set to be either parabolic, 
as expected in low-density environments, or eccentric, 
with an apocenter-to-pericenter ratio of 6:1, as is typi- 
cally seen in cosmological simulations in denser environ- 
ments, i. e., in galaxy group s and in the cores of galaxy 
clusters (|Ghigna et al.lll998ft . 

This paper is organized as follows: In section [2] we de- 
scribe the set-up of the simulated mergers. In section 
[3] we define the tidal parameter and the merger time, 
which are of relevance in the further discussion. In sec- 
tion 14.11 we discuss the strength and lifetime of the tidal 
features, their morphology in section T4.21 and their ex- 
pected B — R color evolution in section [4~3l In section [4~4l 
we estimate the expected relative frequency of elliptical- 
elliptical and elliptical-disk mergers based on statistical 
studies of galaxies in the field. We discuss our results 
and conclude in section [5] 

2. INITIAL CONDITIONS 

Our simulations have been carried out on the high- 
performance Linux clus ter Gonzales using the parallel TV- 
body code PKDGRAV (IStade1 l200llh the gravity modul e 
of the TreeSPH-code GASOLINE (jWadslev et aJ]|200l . 
The galaxy models used in the collisions consist of a 
stellar system embedded in a dark matter halo with an 
NFW (|Navarro et all |1997[ ) profile. The stellar system 
is either a n exponential di sk, a disk with a bulge with 
Hernquist (Hcrnquistl ll990l ) profile, or a stellar model of 
an elliptical galaxy constructed as described b elow. The 
prope rties of the disks are derived analytically (|Mo et al.l 
1998) from the dark matter halo properties. ./V-body re- 
alizations of the specified galaxy models are initialized 
with appropriate velocities and positions of the various 
particles to obtain an equ ilibrium model using a standard 
procedure (Hernquist 1993) widely used in the literature. 
This approach leads to a compound galaxy system that is 
nearly in equilibrium right from the beginning. To min- 
imize the spurious effects due to small departures from 
exact equilibrium we evolved our models for several Gyr 
in isolation before using them in the merger simulations. 
Each of our disk models contains 1-3 million dark matter 
and star particles and the typical force softening is 100- 
300 pc. We explicitly checked that with the chosen mass 
and force resolution our simulations are not significantly 
affect ed by artificial two-body relaxation (|Moore et al.l 
1996), which could cause spurious heating and broaden- 
ing of cold, sharp tidal features. The structural proper- 
ties of our disk models are enumerated in Table [TJ 

Rather than attempting an "ab initio" approach in 
a fully cosmological framework, our massive elliptical 
galaxy is built simply by merging two equal-mass, spheri- 
cal stellar systems embedded in dark matter NFW halos. 
This is the same procedure that has been successfully ap- 
plied to build models of t riaxial dark matter halos start- 
ing from spherical halos ([Moore et al.ll200"4T ) and results 
in a triaxial, boxy stellar system with a small degree 
of rotation. The generated model (and its subsequent 
merger remnants) satisfies the main observational prop- 
erties of luminous ellipticals, including the location on 
the fundamental plane (Fig. [2]) . 

The elliptical galaxy model has a stellar mass of 
2.7 x 10 11 Mq and an effective radius of 9 kpc and we 



set the stellar mass-to-light ratio to 4 in the i?-band, 
thes e numbers being ty pical for massive elliptical galax- 
ies ([Bender et al.l IT992) . The galaxy is resolved with 2 
million dark matter and 2 million star particles, and the 
force softening is 210 pc for dark matter particles and 
60 pc for star particles. Smaller mass elliptical galaxies 
are obtained by rescaling particle masses, positions and 
velocities and softening lengths of the massive elliptical 
based on the general scaling relations for ga l axies formed 
in a cold dark matter model (Mo et al. 1998; Mayer et al. 
l2002h . 

A merger between two initial galaxy models is prepared 
by regarding the galaxies as mass points with masses of 
/i times the galaxy mass. The rescaling factor takes into 
account that in case of close initial positions the galaxies 
are penetrating each other's dark matter halo. The two- 
body problem is solved for the given orbital parameters 
dmit, e, and d per i. Here di n u refers to the initial distance 
between the centers of masses of the galaxies, while d per i 
is the distance at closest approach, the so-called peri- 
center distance. Here e specifies the eccentricity of the 
orbit. The galaxy models are then placed at the location 
of the point masses with their velocities, and additional 
configuration parameters, e.g. the inclination of the disk 
with respect to the orbital plane, are assigned. 

We consider two main classes of orbital configurations. 
The first resembles the nearly-parabolic encounter of two 
galaxies as expected fro m cosmological simulations to oc- 
cur mainly in the field (Kho chfar fc Burkertl l2006). Our 
second class of merger set-ups deals with lower energy, 
more circular orbits such as those expected for galaxies in 
virialized groups and clusters of galaxies (|Ghigna et al.1 
1998). The apocenter-pericenter ratio of about 6 : 1 is 
typical for the latter category of systems. The merger 
parameters are itemized in Table [5J 

3. ANALYSIS 

We define the strength of tidal features by means of a 
tidal parameter 9 based on two-dimensional projections 
of the mass density. We first project the stellar com- 
ponent of our iV-body system onto a two-dimensional 
pixelized grid assuming a constant mass-to-light ratio of 
4 in the i?-band for all stellar particles. This is done 
with an adaptive cloud-in-cell algorithm, where the ker- 
nel scale is chosen according to the local density. The 
images have a resolution of 0.5 kpc pixel -1 , or a scale of 
0'.'27 pixel -1 if we assume that our mergers take place 
at redshift z ~ 0.1. Each image is first corrected for 
(1 + z ) 4 redshift dimming a nd i^-correction (using ZE- 
BRA; Feldman n et al.l [ 206 6), smoothed with a Gaussian 
kernel with FWHM of l'.'l, in order to mimic the see- 
ing conditions, and finally degraded with a background 
noise of 29.5 mag arcsec -2 . This generates our "real" 
galaxy image G. We then mask all pixels that have a 
flux lower than the background noise. Next we apply the 
reverse unsharp masking technique and mask all pixels 
that deviate by more than a factor of 2 after smoothing 
the image with a Gaussian kernel. Pixels neighboring 
a masked pixel are also masked. In the next step we 
fit elliptical isophotcs to our image using the well-known 
ellipse task of IRAF. We fix the center at the posi- 
tion of highest stellar density and allow the ellipticity 
and the position angle to vary freely. This results in the 
model image M. If the two merger constituents have not 
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TABLE 1 

Parameters of the Disk Galaxy Models 



Description 


^200 
(kpc) 


c 


A 


m D 




Za/Rd 


m B 


n/Rd 




N D 


N B 


(kpc) 


ho 
(kpc) 


h B 
(kpc) 


1:4 disk/bulge 


206.3 


12 


0.031 


0.04 


2.59 


0.1 


0.008 


0.2 


1 x 10 6 


1 x 10 5 


1 x 10 5 


0.13 


0.13 


0.13 


1:4 disk 


273 


8 


0.0398 


0.03 


5.64 


0.15 







2.8 x 10 s 


2 x 10 5 





0.3 


0.2 




1:10 disk/bulge 


171.4 


12 


0.031 


0.04 


1.84 


0.1 


0.008 


0.2 


1 x 10 6 


1 x 10 5 


1 x 10 5 


0.13 


0.13 


0.13 


1:10 disk 


200 


8 


0.0398 


0.03 


4.14 


0.15 







2.8 x 10 s 


2 x 10 5 





0.3 


0.2 





Note. — The first column describes the model type and its (stellar) mass ratio with respect to the mass of the massive elliptical used 
in the merger set up. The other columns denote the virial radius, the concentration, the halo spin, the disk mass fraction, the disk scale 
length, the disk height in units of disk scale length, the bulge mass fraction, the bulge scale length in units of disk scale length, the number 
of dark, disk, and bulge particles, and their gravitational softenings. We assume that disk material and dark matter have the same specific 
angular momentum, i.e. jo = mrj. The Hubble parameter is H = 70 km s _ l Mpc -1 . 
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Fig . 1. — Scaling properties of our merger remnants. Blue circles show the properties of observed massive elliptical galaxies (Ben der et all 
1992). The position of our massive elliptical galaxy model is indicated with the green diamond. Other symbols denote our simulated merger 
remnants. The mass ratio (1:1, 1:4, 1:10), the type of the satellite progenitor (elliptical [E], disk+bulge [d/b] or pure disk [disk]) and whether 
the orbit is parabolic (p) or eccentric (e) is indicated for each remnant, (a) Our massive elliptical host galaxy and its remnants obey roughly 
the typical mass-radius relation of m assive elliptical gal axies ( Bender ct al. 1992) . (b-d) The merger remnants lie well within the fundamental 
plane defined by the re-parameters iBender et al.lH992T l rei , K2 and K3. They span an orthogonal coordinate system in the three-space of 
the logarithm of the central velocity dispersion, the logarithm of the effective radius, and the mean surface brightness within the effective 
radius. Here Ki, K2, and re3 are approximately proportional to the logarithm of the stellar mass, the mean surface brightness, and the 
logarithm of the mass-to-light ratio, respectively. 
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TABLE 2 

The Orbital Parameters of the Various Merger Setups 







r init 


r peri 






Description 


Orbit 


(kpc) 


(kpc) 


e 


H 


(1:1) E + E 


P 


700 


50 


0.95 


1 


(1:4) E + E 


P 


518 


18 


0.99 


1 


F! -i- F, 
^r.±uy 1 i ^ Hi 


P 


518 


18 


0.99 


j 


(1:4) E + disk/bulge, 












rpcri 9 kpc 


P 


518 


9 


0.99 


1 


(1:4) E + disk/bulge 


P 


518 


18 


0.99 


1 


(1:4) E + disk/bulge, 












rperi 36 kpc 


P 


518 


36 


0.99 


1 


(1:4) E + disk 


P 


518 


18 


0.99 


1 


(1:10) E + disk/bulge 


P 


518 


18 


0.99 


1 


(1:10) E + disk 


P 


518 


18 


0.99 


1 


(1:1) E + E 


c 


250 


35.7 


0.75 


1 


(1:4) E + E 


c 


150 


15 


0.818 


0,1 


(1:10) E + E 


c 


150 


15 


0.818 


0.4 


(1:10) E + disk 


c 


150 


15 


0.818 


0.4 



Note. — The columns contain: the identification of the respec- 
tive merger, whether the orbit is parabolic (p) or eccentric (e), the 
initial distance, the Keplerian pericentric distance, the Keplerian 
eccentricity, and the mass reduction factor fi of the major progeni- 
tor (se e text). The first simulation is taken from Kazantzid is et al.l 
(2005). The inclination angles between orbital and spin plane 
for the major and minor participant are chosen randomly in our 
simulations but fixed for a given orbital type (parabolic or eccen- 
tric) with the constraint that coplanar or polar configurations are 
avoided. The pericentric arguments are also chosen randomly. 

merged we apply the described technique to each object 
individually and add the model images to create M . A 
distortion image F is created by first dividing G by M 
pixelwise and then convolving it with a quadratic median 
filter of 5 x 5 pixels to reduce pixel-to-pixel variation. The 
tidal parameter 9 is now defined as 9 = (\F— 1|). It mea- 
sures the (absolute) deviation of the fractional residuals 
over the nonmasked pixels of the image. 

As an operational definition of weak and strong tidal 
features we use values of 9 = 0.10 and 0.15 as our "weak" 
and "strong" tidal feature limits. The former limit cor- 
responds to weak, barely visibly tidal features and the 
latter to strong, easily identifiable tidal features. For 
comparison, we measured a value of 9 = 0.085 ± 0.002 
for our massive elliptical galaxy model before the inter- 
action begins. The lifetime of weak/strong tidal features 
is understood as the characteristic time during which the 
median of 9 over the various projections is larger than the 
weak/strong feature limit. The result of this analysis is 
shown in Table [3] 

The morphological properties of merger progenitors 
and remnants are analyzed from the projected images. 
Using the ellipse task of IRAF the cllipticity e and the 
isophotal shape 100a4/a are determined as a function of 
semi-major axis length a. The effective radius is esti- 
mated as a product between semi-major axis length ah 
containing half the light and (1 — e/j) 1 / 2 , where eh is the 
cllipticity measured at ah- The central velocity disper- 
sion is measured within a circular aperture of lkpc radius 
mimi cking observations (jBender et al.lfl992l : lDavies et al.l 
1987). In Tabled we list the median structural and kine- 
matic properties averaged in the range of 0.8 - 1.2 times 
the half-light semi-major axis as obtained from at least 
50, but typically more than 200 random, projections. 

We define the merger time based on the particle data. 



In particular, we keep track of the distance between the 
star particles with the highest density of each merger 
constituent and define the time of merging as the time 
when this distance drops below 5 kpc (Fig. [3]). 

The overall merging time is shorter for more massive 
companion galaxies. It is determined by a combination 
of the dynamical frictio n time with the tim escale of mass 
loss via tidal stripping (|Taffoni et al.ll2003T ) . The dynam- 
ical friction timescale is longer for a lighter companion 
and t idal stripping slow s it down further (|Colpi et al.l 
Il999t iTaffoni et aLll2003l) . More than half of the mass 
of the dark halo is typically stripped along the first or- 
bit, while the stellar component loses mass more gradu- 
ally since it sits deeper in the potential well of the halo. 
When measured from first pericenter the merger takes 
on parabolic orbits < 1 Gyr for equal mass galaxies, 
and between 7 and up to 12 Gyr for a 1:10 mass ra- 
tio. The merging time is generally faster on eccentric 
orbits because they start with a smaller orbital time in 
our initialization, which yields a shorter dynamical fric- 
tion time. Finally, pure disk companions are disrupted 
much faster than disks with bulges or elliptical compan- 
ions due to their lower central density, which results in 
a smaller tidal radius and a m ore impulsive respon se to 
the tidal shocks at pericenters (iGnedin et al.lfl999T ). 

An observer would define a merger as complete when 
the individual merger components could not be distin- 
guished anymore. Using our mock observational setup 
we find that equal mass mergers happen on such a fast 
timescale that the particle-based merging time and the 
observational merging time coincide. A similar statement 
holds for the simulations involving only elliptical galaxies 
since in this case the core remains extremely compact and 
reaches the inner kiloparsec of the host galaxy without 
being destroyed. For our 1:10 merger involving a pure 
disk on the other hand there is a significant difference 
of ~ 200 Myr; i.e. the progenitors have not yet merged 
according to the particle definition, but the satellite is so 
strongly disrupted that it is difficult to recognize it as a 
single individual object. 

4. RESULTS 

4.1. Strength and lifetime of tidal features 

Fig. 14.11 shows an example of the dynamical evolution 
of one of the elliptical-disk merger simulations. The mass 
ratio between the satellite disk and the main elliptical 
is 0.1. Stars belonging to the disk of the small satel- 
lite are stretched in tidal streams at the first pericentric 
passage. These tidal features, which originate from the 
companion's disk component, survive for a few Gyr after 
the merger, in contrast with the short-lived debris that 
is produced in elliptical-elliptical mergers. We quantify 
the strength of the tidal features by means of the tidal 
parameter 9; this measures the light excess in the tidal 
tails relative to a smooth elliptical-isophote model of the 
galaxy, see section [3] 

Mergers between elliptical galaxies, bona fide dry 
mergers, produce tidal features that are short lived; 
in particular equal-mass mergers do not show strong 
features for more than 100 Myr (Table [3]). Mergers 
with disk-dominated satellites, instead, produce not only 
stronger but also longer lived features. Strong, easily 
identifiable features (9 > 0.15) last for more than 1 Gyr, 
while weak, barely visible features (9 « 0.1) persist for 
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TABLE 3 

The Lifetimes in Strong, Intermediate and Weak Tidal Features for the Different Simulations 



Description 


Orbit 


Weak (0 > 0.1) 
(Mvr) 


Intermediate (6 > 0.13) 
fMvr) 


Strong (9 > 0.15) 
fMvr~) 


1:1 E + E 
1:4 E + E 
1:10 E + E 


P 
P 
P 


1359+ 2 ? | 

1339_32 2 


74.76+H 2 

485.3_«*;_ 


2s.oiiJ2:I 

91 R 9+29.8 
z±o.z_ g 32 

0+° 


1:4 E + disk/bulge, rperi 9 kpc 

1:4 E + disk/bulge 

1:4 E + disk/bulge, rperi 36 kpc 

1:4 E + disk 

1:10 E + disk/bulge 

1:10 E + disk 


P 
P 
1) 
P 
1) 
P 


3150_2„ 
2692l|3.2 

> 1980 
3196+1?;? 
413.2l| 5 5 « 
2028l 7 3;i 


1472+ 82 - 4 
2109+?|? 

lesel^ 3 

1906+21 1 
0+0 

1007+«-j 


1508l 6 3 ^ 
1287±p;| 
1537+ 37 ; 8 

0+° 
610.5+^8-j 



1:1 E + 
1:4 E + 
1:10 E 


- E 
E 
+ E 


c 
e 
c 


28.04±f 7 7 8 
12191S.! 
916.8±^ 


+0 
U -0 
441.6t 3 2 6 4 l 

6.087i^ 9 


0+0 

U -0 
o 97 +43.3 
oz -49.3 

0+0 


1:10 E 


+ disk 


e 


343112.9 


1344± 6 32 : 2 8 


650.31^:2 



Note. — The first column denotes the mass ratio and the type of the merger progenitors. The single letter in the second column 
specifies whether the binary system is initially bound and the orbit is eccentric (e) or the orbital energy is zero and the orbit is initially 
parabolic (p). The error limits correspond to lower and upper qua rtiles derived by adding random values with mean zero and a spread 
determined from the lower and upper quartiles shown in Figure |4~T1 to the original tidal parameter curve and measuring the tidal-lifetime 
for each randomization. The values of corresponding to weak and strong features are chosen in a way such that weak features are barely 
noticeable under visual inspection of the mock image while strong features are outstanding emergences of tidal distortions. 



TABLE 4 

The Kinematic and Structural Properties of the Massive Elliptical and the Merger Remnants 



Description 


Orbit 


M e 

(1O 1O M ) 


(kp°c) 


(km s" 1 ) 


e 


100a 4 /a 


Massive elliptical (E) 




13.84 


q 24+0-23 


248_ 9 4 8 


n noo + 0.039 

U.OOO_ 98 


— 2.58l .47 


1:1 E H 
1:4 E H 
1:10 E 


- E 

- E 
+ E 


P 
P 
P 


27.69 
17.3 
15.23 


16.8+ij 
13.9t°/f 

ii O+0.49 
11.S_0.37 


296+JI 
277l?g 
283_ 6 i ; 7 


4 7+0.079 
u -^'-0.16 
Q 4 1 +0.061 
U.CS41_o 077 

261 + 043 
U.ZDl_o 057 


-2.4llJ.46 
0.585±°i| 
_n q 7 +0 34 


1:4 E + disk/bulge, rperi 9kpc 
1:4 E + disk/bulge 
1:4 E + disk/bulge, rperi 36 kpc 
1:4 E + disk 
1:10 E + disk/bulge 
1:10 E + disk 


P 
P 
P 
P 
P 
P 


17.7 
17.7 
17.7 
17.3 
15.22 
15.21 


10 9+0.61 
lo ' z -0.79 
in O+0.79 
lz -°-0.82 

12 Q+ - 67 
J^-y_0.95 

14.O_0.88 

11 1+0.57 

^^-o.ss 

ii fi +0.59 
li ' D -0.55 


266+11 
26611 9 
267+?2 
2641-8 
2771 It 

2™±"s 


n o 9 c+0.053 
U.JZO_ 086 

S69+°- ()G1 
U.JDy_ 083 
oqi +0.057 
U.O»1_ 088 
n qoo+0.06 
U.OOJ_o 086 

U.JUO_ 091 
9 4c; +0.056 
U.Z4O_ 062 


o.i33l°;S 

0.454l0.3 2 4 
fi4 +0.43 
U - D *-0.31 
Q9c;+0.38 
u.ozo_ 29 

0.789ir 4 3 8 
-0.854l°;35 


1:1 E H 
1:4 E H 
1:10 E 


- E 

- E 

f E 


c 
c 
c 


27.69 
17.3 
15.23 


14.8+1- 8 

ii fi +0.91 
11.D_0.61 
1n O+0.73 
iU - S -0.66 


2891^2 
287l2l 
2771^ 


0.34ll°13 

n 9HS+0.082 
0.z»8_o 064 
n 9^7+0.064 
U.ZO(_ Q 06 


0.0515tg-2? 
0.7591^5 
-0.7571°;!? 


1:10 E + 


disk 


c 


15.21 


10 fi+ - 75 
1U.D_ 59 


270l 9 2 3 


24+ - 06 
u - z ^-0.082 


-0.277lg'.2? 



Note. — The first three columns denote the merger progenitor or remnant, respectively, the orbital setup (parabolic or eccentric), and 
half of the total system's mass. The values of the following columns express the median values over random projections. In particular, the 
effective radius, the central velocity dispersion within 1 kpc, the ellipticity at effective radius, and the isophotal shape are given. Error 
limits refer to lower and upper quartiles. The properties of the mass-rescaled versions of the elliptical progenitor (E) can be obtained from 
the scaling relations for the coordinates and velocities. 
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Fig. 2. — Distance of the star particles with highest local density from each progenitor as function of time. The first time that the 
latter distance drops below 5kpc and (except for noise) remains below this threshold defines the merging time. Note that mergers of a 
given mass ratio between their stellar components but with different satellite type may have different total mass ratios (the total mass 
including dark matter), (a) Mergers on parabolic orbits between a massive elliptical and pure disk satellites (disk) or satellites with disk 
and bulge (d/b). (b) Mergers between elliptical galaxies (E) on parabolic orbits, (c) Mergers between elliptical galaxies on eccentric orbits. 
(d) Merger between a massive elliptical and a 10 times less massive disk satellite on an eccentric orbit. 
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Fig. 3. — Evolution of merging galaxies. Surface brightness maps showing the stellar component at different stages of the simulation of 
a merger between a massive elliptical galaxy and a bulgeless disk with 1/10 of the mass of the accreting elliptical; the systems are on an 
eccentric orbit. Top: Merging systems 200 Myr after the beginning of the simulation (left) and 800 Myr later on their way to the second 
pericentric passage (right). At this point the disk is already heavily damaged and tidal streams appear. Bottom: Remnant 500 Myr (left) 
and 1 Gyr (right) after the merger. The visible tidal streams stem from the disk and begin to develop at the first pericentric passage. They 
survive up to a few Gyr after the merger. All images have a physical size of 300x300 kpc. The surface brightness ranges from 28.5 mag 
arcsec -2 (black) to 18 mag arcsec -2 (white). Noise is not added to these images. 



up to 3 Gyr. We explain this difference with the smaller 
velocity dispersion of the disk material, which makes up 
for most of the stellar tidal debris even in the galaxies 
with a bulge component, which implies that phase mix- 
ing with the background mass distribution occurs on a 
longer timescale and that tidal features are thus longer- 
lived compared to dry mergers. 

We find that the time dependence of the tidal param- 
eter can be well fit with a declining exponential of the 
form 6(t) = Ae-t^ + B fFig.P]). The additive constant 
B is kept fixed at B — 0.085, which is the value that we 
measure for our elliptical galaxy model in isolation and 
only A and r are fitted. The product of A and r, which 
represents the area under the curve of the tidal parame- 
ter once the offset B is removed and is hence a measure 
of the strength of the tidal features, discriminates well 
between dry mergers and mergers involving disks. Each 
of our dry merger simulations shows At < 0.07 Gyr while 
all our disk simulations have At > 0.1 Gyr, with the no- 
table exception of the 1:10 merger involving a disk+bulge 
system. 

The strength of tidal features originating from a merger 
with a pure disk satellite and from a merger with a 
disk plus bulge satellite are found to be of comparable 



strength for a mass ratio of 1:4 but very different for 
a mass ratio of 1:10. In particular, the merger with a 
disk plus bulge system does not show strong tidal fea- 
tures at merging time as opposed to the strong and long- 
lived features in the 1:10 merger with a pure disk satel- 
lite. We explain this difference as follows. In our sim- 
ulation the bulge-less satellite is disrupted before prop- 
erly sinking to the center of the giant elliptical. Hence, 
the merging time, which is determined by the minimum 
of disruption time and tidal friction time, is set by the 
smaller of the two, which in this case is the tidal dis- 
ruption time. On the other hand galaxies containing a 
bulge are more concentrated and hence less susceptible 
to both tidal stripping and tidal heating. This in effect 
increases the tidal disruption time and consequently the 
merging time, which is no w essentially determ ined by the 
dynamical friction time (Taffo nTet al.ll2003T ). The con- 
centrated, cold tidal features in our simulations originate 
from the stellar disk and not from the kinematically hot 
bulge stars. The features are mainly created at peri- 
centric passages and at the final coalescence of the disk 
component with the elliptical galaxy (as tidal forces scale 
with the cube of the inverse distance). Hence, the pres- 
ence of a bulge introduces a significant time span between 
the creation of the tidal features and the final merger of 
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time since merger in Gyrs 

Fig. 4. — Tidal parameter 9 as function of time elapsed after the merger. For each merger set-up and time step, the tidal parameter is 
calculated over 20 random projections. The colored dashed lines show its median values, and the shaded area defines the lower and upper 
quartile. Colored dotted lines indicate exponential fits to the declining tidal parameter, (a) Parabolic encounters involving disk galaxies. 
Shown are a 1:4 merger with a pure disk satellite (red), a 1:4 merger with disk+bulge satellite (black), a 1:10 merger with a pure disk 
satellite (blue), and a 1:10 merger with a disk+bulge satellite (green). The strongest features (9 > 0.15) are visible for the most massive 
disks and, in this example, last for about 1.5 Gyr. (6) Dry mergers between elliptical galaxies on parabolic orbits with stellar mass ratios 
1:1 (black), 1:4 (red) and 1:10 (magenta). Strong features are visible for only a very short time, (c) Encounters on eccentric orbits. Black, 
red, and magenta curves denote elliptical-elliptical mergers as in (6). In addition, a 1:10 merger with a disk galaxy is shown (blue). 



the dense core of the satellite with the center of the giant 
elliptical. The produced tidal features may thus become 
severely dimmed or even invisible at merging time. In 
the case of the 1:4 mass ratio this time span is small for 
both the satellites with bulge and without because for 
companions of such a high mass the dynamical friction 
timescale is always dominant in setting the merging time 
(the dynamical friction time is only mildly sensitive to 
the internal structure of the secondary galaxy) . 

We also study the dependence of the tidal debris on the 
orbital parameters of the encounter for a given structure 
and mass ratio of the two galaxies. Rerunning the 1:4 
merger involving a disk plus bulge satellite on a parabolic 
orbit with twice and with half the Keplerian pericenter 
of 18 kpc does not change the tidal strength and life- 
times significantly. For eccentric bound orbits, which 
is the more natural assumption if the merger occurs in 
a virialized group rather than in the field, tidal debris 
lasts slightly longer than for the corresponding mergers 
on parabolic orbits. The reason is that disruption occurs 
deeper in the central part of the potential of the primary 
for such eccentric orbits. Closer to the center of our 
elliptical galaxy model we find the inertia tensor to be 
less triaxial and the effect of orbital precession, which 
wraps ar ound and phase m ixes tidal streams in triaxial 
systems (jMaver et al.l 120 02). is reduced. 

4.2. Morphology of the tidal debris 

A striking difference between elliptical-elliptical and 
elliptical-disk mergers is manifested in the morphology 
of the resulting tidal tails (Fig. [5}. Mergers between 
two kinematically hot early-type systems produce dif- 
fuse plumes of stars, which vanish rapidly (|van Dokkuml 
[200alBell et"aLl[2006h . In some of our simulated mergers 
between early-type galaxies of unequal mass we observe 
in addition straight, elongated tails, which also vanish 



rapidly. On the other hand, we find that mergers that in- 
volve a dynamically cold disk are able produce strong and 
long-lived tidal features, with loops and ar ms resembling 
those observed around bright ellipticals (|van Dokkuml 
f2005h . 

In order to estimate the fraction of systems with broad 
tidal features versus features with loops and arms we vi- 
sually classify the morphology of the tidal debris in the 
sam ple of 86 nearby b ulge-dominated red field galax- 
ies (Ivan Dokkuml 120051) . Twenty-five of the selected 
galaxies are classified as possessing no observable tidal 
features, 27 objects are classified as showing weak fea- 
tures, 14 objects show strong features, and 20 objects 
are currently undergoing a merger. For six galaxy pairs 
the merger companions are included in the sample of the 
20 merging galaxies, leaving 17 unique red merger pairs. 
We find that of 57 unique sample galaxies for which signs 
of interactions are detected (out of the total sample of 
83=86-3 bulge-dominated, red field galaxies), only nine 
systems show the "broad fans" that are expected from 
a very recent elliptical-elliptical merger. In contrast, 19 
systems clearly show tidal arms and loops. This statis- 
tics should of course be taken with care, in the light of 
the 29 cases in which the tidal features are too faint for 
a decisive classification of their morphology. Hence, our 
qualitative visual inspection confirms that the majority 
of the observed tidal features in elliptical galaxies are 
indeed consistent with mergers involving disk satellites 
rather than bona fide dry mergers. 

4.3. The color of the tidal features and the merger 
remnant 

In the following discussion we assume Johnson Bj 
and Cousins Rq filter bands. Magnitudes and colors 
are normalized to Vega. Our simulations show that 
the tidal arms and loops stem mainly from the regions 
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Fig. 5. — Morphologi cal sig natures of mergers, (a) Observed image of an elliptical galaxy that shows red tidal tails at large galactocentric 
distances (van Dokkum 2005). (b) Mock image of a simulated 1:4 merger between an elliptical galaxy and a disk+bulge system; the snapshot 
is taken 600 Myrs after the merger, (c) Mock image of a simulated 1:10 merger between an elliptical galaxy and a disk; the snapshot is 
taken 500 Myrs after the merger, (d) Mock image of a simulated 1:4 merger between elliptical galaxies, 380 Myrs after the merger. There is 
a striking difference between the tidal features originating in mergers between ellipticals and those involving a disk-dominated companion, 
independent of the mass ratio and orbit characteristics. While all mergers can lead to shells and diffuse features, only the mergers involving 
disks show strong tidal arms and loops, similar to the observed features around bright ellipticals (see [a]). Panel a has a size of 2.'5x2f5, 
corresponding to a physical extension of about 280x280 kpc for an object at redshift z ~ 0.1. The mock images (b-d) have a physical 
dimension of 300x300 kpc and hav e been generated with resolution and noise properties so as to match the observational data used for 
the comparison ( van Dokkum 20055). The bright internal regions of the images are saturated, in order to increase the contrast on the tidal 
features. 



of the companion's disk at about 2 — 5 radial scale 
lengths from the center (i.e., at about the solar neigh- 
borhood in the Milky Way). The observed (B - R) 
color of M31 at such distances, for example, is about 1 .4 
mag (fWalterbos fc KennicuttJll988HTamm et al.ll2007t ). 4 
This is only ~ 0.1 mag bluer than the typical color of gi- 
ant ellipticals at 30-50 kpc from the center, where the 
tidal tails are observed. Indeed, the cores of nearby 
gian t ellipticals ha ve typically (B — R) ~ 1.6 mag, 
(e.g. lMichardlll999h . and, given the typic al color gradi- 
ents t hat are observed for giant ellipticals (|Peletier et al.1 
I1990I) . the color at 30-50 kpc is typically (B - R) ~ 1.5 
mag (Fig. [6|). 

To stress this point further, we compile in Fig. [TJi the 
(B — R) colors as a function of sc a le len gth from the 
sample of Ide Jong fc van der Kruitl (|1994| ); the sample 
contains 86 face-on nearby disk-dominated galaxies. We 
note that (a) the typical (B — R) color of disk galaxies 

4 The stellar mass ratio between M31 and our elliptical galaxy 
model is about 1:4. Less massive disks will be on average less 
luminous, more gas-rich, and bluer. 



between 2 and 5 radial scale lengths is ~ 1.2 mag; and 
(b) the 1 a fluctuations for individual galaxies are of the 
order of 0.25 mag. For a direct comparison, we show in 
Fig. EJd the (B — R) color of the sampl e of 39 giant el- 
liptical galaxies of F eletier et al.l (|1990f ). At about two 
half-light radii, i.e., at the distance at which tidal loops 
and arms are observed, the typical (B — R) color of el- 
lipticals can be seen to be ~ 1.5 ± 0.1 mag. Hence, the 
(B — R) color difference between the disk stellar popu- 
lation responsible for the tidal arm and the background 
population of an elliptical galaxy is of the order of 0.3 
mag, and its typical scatter is of the same order. This 
difference in colors of a few tenths of a magnitude is 
already quite small but has to be viewed as a conser- 
vative iminbCT^^n_Jact J _the_diskgalaxies in the sam- 
ple of Ide Jong fc van der Kruitl (fl99l are normal star- 
forming disk galaxies, while in our case the disk satellite 
will have the star formation quenched in the outer region 
as soon as it completes the first close approach with the 
primary elliptical as a result of gas stripping. As seen in 
hydrodynamical simulations tides will remove the gas at 
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Fig. 6. — Dependence of the (B — R) color on radius mea- 
sured in physical kiloparsecs. Lines denotes the average (B — R) 
color in the radial bin of the r espective sample of spiral galaxies 
(dc Jong & van dcr Kruit 1994, blue solid line) and giant elliptical 
galaxies llPeletier et alj|199& red dashed line). Error bars denote 
the 1 cr error of the average (B — R) color. Circles indicate the 
median (B — R) color in the respective radial bin. The shaded area 
shows the 1 a scatter for individual color measurements. 



least down to the same radius at which stars are removed 
(which is down to 2.5 disk scale lengths in our simula- 
tions) and ram pressure in the gaseous halo that the giant 
elliptical will presumably posses s will enhance stripping 
even further ([Mayer et al.l l2006.. Tidal tails originating 
from 2-5 disk scale lengths would thus be indeed red - 
not blue - features. 

This conclusion is strengthened by an analysis 
con ducted with stell ar population synthesis mod- 
els (|Bruzual fe Charlotf2 003 ) . In Fig. M we plot the B-R 
color as a function of epoch for simple stellar populations 
(red solid line) . For our elliptical of mass ~3x 10 11 
we assume overall solar metallicity (|Kodama et all l998). 
We note that the (B — R) color varies between -0.2 and 
1.6 mag, i.e., it spans more than 1.8 mag, when moving 
from very young (< 10 7 yr) to very old (12 Gyr) stellar 
populations. We model the (B — R) color of a ~ 3 x 10 10 
Mq disk galaxy by assuming an exponentially declining 
star- formation rate with a 4 Gyr e-folding star fo rma- 
tion timescale ([Lilly et al.lll998llBell fc de Jondl2000l :blue 
solid line). We find that the (B — R) color of our disk is 
typically only ~ 0.4 mag smaller then the (B — R) color 
of the bulk stars of the massive elliptical. We also con- 
sider a set of models with e-folding times ranging from 
1.5 Gyr to infinity (i.e. constant star formation rate) and 
metallicities varying from half-solar to solar. These mod- 
els predict a spread in (B — R) color of about 0.4 mag. 
The horizontal dashed lines indicate the lower and upper 
quartile of the (B — R) color of face-on disks at about 3 
disk scale lengths (Fig. [7^). Models in which the star for- 
mation in the outer disk is suppressed as a result of gas 
stripping are also included (blue dashed lines). The con- 
sidered timescale of quenching of star formation previous 
to observation of the tails ranges from 4 Gyr (upper blue 
dashed line) to 1 Gyr (lower blue dashed line). These 
models exhibit a (B — R) color between 1.26 and 1.38 
mag, i.e, only 0.1 - 0.25 mag less than the color of the 



bulk stars of the elliptical at the radius at which tidal 
tails are observed. Note that 1 Gyr is a conservative 
timescale, since material from the outer disk is stripped 
well before the merger is completed and the strong tails 
remain visible for up to 1.5 Gyr in disk-elliptical mergers. 
Our stellar sy nthesis models employ a Chabrier initial 
mass function (|Chabrieril2003l ). 

We therefore conclu de that in the qualitative analysis 
of Ivan Dokkuml (|2005[ ) the colors of tidal debris originat- 
ing from the disruption of a disk companion would be 
within at most a few tenths of magnitude from the col- 
ors of the background stars of the host elliptical galaxy 
at those galactocentric distances and would thus be iden- 
tified as "red" features. 

We note that in principle on-going star formation dur- 
ing and after the completion of the merger could lead 
to tidal features with bluer colors. This would weaken 
our conclusion that disk-elliptical minor mergers are re- 
sponsible for red tidal features. However, some observa- 
tions suggest that blue tails should exist only for a very 
short time after the merger and should redden quickl y 
(|Schombert et all Il990t iHibbard fc van Gorkoml 1 19961) . 
Nevertheless, a further study of this issue is necessary 
and will require the use of hydrodynamical simulations 
with a realistic model of star formation. 

In Fig. [9] we investigate the color of the massive ellip- 
tical after a 1:10 merger with a disk galaxy that contains 
10% of its mass in the form of gas. Within less than 300 
Myr the B-R colors of the re mnant are red enoug h to fall 
into the selection window of Ivan Dokkuml (j2005f ). 

4.4. The fraction of tidal features originating from 
mergers with disks 

For any given mass ratio, we calculate the average rela- 
tive lifetimes ft of the weak (Tabled tidal features gener- 
ated respectively by elliptical-disk and elliptical-elliptical 
mergers (averaged over various simulation set-ups). This 
represents a conservative lower limit on the relative life- 
time, since stronger tidal features would have even larger 
relative lifetimes. We estimate that the average disk- 
elliptical-merger lifetime is a factor ~ 2.3 larger than the 
elliptical-elliptical- merger lifetime for a 1:4 mass ratio; 
this factor becomes > 4 for a 1:10 merger. 

We use t he i3-band luminosity functions 

of iDriver et al.l (|2006l ) to give a rough estimate of 
the relative mass fraction of early-type versus disk 
galaxies (see Fig. [T0k.. Those authors classify by visual 
inspection galaxies into E/S0, Sa/b/c, or Sd/Irr/Pec. In 
the following analysis we identify the first class (E/S0) 
with early-type galaxies and the latter classes with 
disks. Our results do not depend on whether we include 
or exclude the Sd/Irr/Pec class because within our 
considered mass range such galaxies do not contribute 
significantly to the mass function. In order to obtain 
estimates for the mass functions near ~ 5 x 10 10 we 
divide the provided luminosity functions by a _?-band 
mass-to-light ratio of 4.5 (E/S0), 2 (Sa/b/c), and 1 
(Sd/Irr/Pec). These mass-to- light ratios are consistent 
with the considered range o f ( B — R) colors deriv ed from 
stellar synthesis modelling ; Bell fc de JondfeOOll ). 

Our elliptical galaxy host has a stellar mass of 2.7x 10 11 
Mq, which corresponds to 2.3M* [with M * ~ 1.2 x 10 11 
Mq pr ovided by the Schechter function fit of lDriver et al.l 
(2006) to the luminosity function of early-type galaxies 
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Fig. 7 . — Dependence of the (B — R ) color on radius measured in units of galaxy's scale lengths, (a) The (B — R) color of spiral 
galaxies (dc Jong & van dcr Kruit 1994) as function of radial disc scale length. Disc scale lengths are derived by a linear fit to the 
surfa ce brightness curve between 10 and 50 arcsec. (b) The (B — R) color of giant elliptical galaxies as function of half-light radius (data 
from lPeletier et al.|[T990h . The major axes provided in the original reference have been transformed to half-light radii using the ellipticity 
at 10 arcsec. Symbols and shaded areas are as in Fig. [6] 



and assuming (M/L)b — 5.5]. For a mass ratio of 1:4 and 
1:10, respectively, we obtain a relative number density of 
disk versus early-type galaxies of 0.60 and 1.45 (see inset 
of Fig. [TUb ) . We have applied the same analysis to the 
mass functions derived by other authors who split galax- 
ies acco rding to di f ferent criteria. From the mass func- 
tions of iBell et al.l (|2003f ) we compute relative number 
densities of 0.64 and 1.32 for 1:4 and 1:10 mergers, re- 
specti vely. The mass functions provided in lBaldrv et al.l 
(2004) lead to relative number densities of 0.81 and 1.05. 

We shall now assume that the relative merger fre- 
quency is proportional to the relative number density 
of disk versus early-type systems. The relative fraction 
of tidal features will thus be given by the product of 
the relative lifetimes of the tidal features with the rel- 
ative number densities of such galaxies. Hence, for a 
2.7 x 10 11 Mq elliptical galaxy the relative fraction of 
tidal tails generated by encounters with disks, relative to 
encounters with other early types, is found to be > 1.4 
for a 1:4 mass ratio, and > 4 for a 1:10 mass ratio. There- 
fore, a lower limit is that more than 60% (mass ratio 1:4) 
or 80% (mass ratio 1:10) of the visible tidal features will 
stem from a merger with a disk-dominated satellite. 

For a fixed mass ratio, decreasing the mass of the main 
elliptical galaxy (i) enhances relative lifetimes, since only 
the strongest features remain visible; and (ii) increases 
the fraction of blue versus red galaxies. Both effects tend 
to increase the relative fraction of tidal features. We can 
estimate the relative fraction f n {R) of tidal features after 
integrating over all possible masses of the elliptical host 
galaxy for a fixed mass ratio R < 1. We will assume that 
the relative lifetime ft{R) only depends on R and not 
on the mass itself. The relative fraction of tidal features 
f(R) is, in this simple approximation, given by: f(R) = 
f n (R)ft(R). In particular, 



f(R) 



MR) 



n D (mR) 
riE\Tn) ; — —dm 



f tie (m) dm J ay ' nE(mR) 
where tie and tie denote the number density of early- 
type and disk galaxies, respectively. The integrals are 



well-defined and need only be evaluated near the charac- 
teristic mass M* of the early-type galaxy population (e.g. 
the integration between 10 9 and 10 12 Mq is accurate up 
to a few percent for all relevant mass ratios). Fig. [TUb 
sh ows fn(R) c a lculat ed from the luminosity functions 
of iDriver et alj (|2006f ) after converting them into mass 
functions as described above. We obtain /„(1 : 4) = 2.03 
and /„(1 : 10) = 2.96. Hence by integrating over all 
plausible masses of the elliptical galaxy host we estimate 
that the relative tidal fractions are > 2.8 for a mass ratio 
of 1:4 and > 12 for a mass ratio of 1:10. 

5. DISCUSSION AND CONCLUSION 

We have found major differences between the elliptical- 
elliptical and elliptical-disk mergers in both morphology 
of tidal tails and their timescales of observability. Only 
mergers that involve a dynamically cold disk produce 
strong tidal features with loops and ar ms resembling 
those observed around bright ellipticals (jvan Dokkuml 
2005); these features last much longer (1-2 Gyr) than 
the smooth features produced in hot-hot mergers. 

Dynamical friction and tidal stripping determine the 
specifics of the merger evolution. For a given satellite 
galaxy type, either disk or elliptical, the tidal tails gen- 
erated by less massive satellites extend to larger distances 
from the center of the accreting elliptical, as the smaller 
the satellite, the longer it needs to sink to the bottom 
of the gravitational potential. However, the smaller the 
satellite mass, the weaker the tidal debris, since less mass 
is available to produce the tidal features. We find that 
satellite masses of order one fourth of the mass of the par- 
ent elliptical exhibit the most long-lived tidal features. 

Since mergers with disk-dominated, relatively low mass 
companions produce the strongest and longest-lasting 
tidal streams, such mergers are expected to be the typical 
events that are revealed in observations of real galaxies. 
In first approximation, the probability that the main el- 
liptical merges with a companion of a given type will be 
proportional to the number density of galaxies of that 
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Fig. 8. — Evolution of (B—R) co lor in synthetic stellar population 
models (Bruzual & Chariot 2003). The (B — R) color as function 
of age for a simple stellar population of solar metallicity (red solid 
line). The blue solid line shows the (B — R) color of a composite 
stellar population with exponential declining star formation rate, 
a 4 Gyr e-folding time, and half-solar metallicity. The shaded 
area brackets the (B — R) colors once we vary the e-folding time 
between 1.5 Gyr and infinity, i.e. a constant star formation rate, 
and the metallicity between half-solar and solar. The colors are 
assumed to be measured 12 Gyr after the initiation of the star 
formation (vertical dot-dashed line). The horizontal dashed lines 
indicate the lower and upper quartiles in (B — R) color at 3 disc 
scale lengths from the center of the disk (from Fig. [7ji). Dashed 
blue lines indicate the (B — R) color of a stellar population with 
initially constant star formation that is then quenched 4, 2 or 1 
Gyr before the observation. The inset shows the (B — R) color 
evolution restricted to the last 6 Gyr previous to observation. O ur 
models assume a Chabrier initial mass function (Chabrier 2 01331) . 

type. The ratio of observable tidal features that are ex- 
pected from mergers involving a disk galaxy, relative to 
mergers that involve another elliptical, can thus be ap- 
proximated by the product of the ratio between the re- 
spective tidal lifetimes and the ratio of the respective 
number densities. Using publish ed luminosity fun ctions 
for early-type and disk gal axies dDriver et 3.1.1120061) and 
typical mass-to-light ratios (|Bell fc de Jong||200lD . we es- 
timate that, at the mass scale of the host galaxy and for 
a 1:4 merger, > 60% of the observed tidal features stem 
from an encounter with a disk-dominated system; the 
fraction increases up to > 80% for 1:10 mergers. We 
thus expect that the merger fraction that is estimated 
from observed tidal features around bright ellipticals re- 
flects primarily the fraction of spheroids that underwent 
(relatively minor) accretion events of dynamically cold, 
disk galaxies over the past several billion years, rather 
than very recent major mergers between elliptical galax- 
ies. 

Our simulations show that the observable tidal debris 
mostly originates from regions of the disk corresponding 
to about 2-5 exponential scale lengths. Using observed 
photo metric profiles of spiral (Ide Jong fc van der Kruitl 
Il994h and elliptical (iPeletier et al.l I1990D galaxies and 
stellar synthesis models ()Bruzual fc CharlotJ 120031 ). we 
find that the B — R color of the tidal features is rather 
similar to the color of the main accreting elliptical at sev- 
eral half-light radii from its center, where the tidal tails 
become visible. We estimate that both the average B — R 
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Fig. 9. — Color-magnitude selection of l umino us red galax- 
ies. The selection window of /van Dokkum (2005) is indicated 
by dott ed lines. The dashe d black curve corresponds to a 3L* 
galaxy (Blanton et al. 2001). The solid lines show the expected 
colors and magnitudes as a function of time (0.1-9 Gyr) and for 
different redshifts (from 0.04 to 0.2) of an old stellar population 
of mass 2.7 X 10 Mq after merging with a 10 times less massive 
disk. This disk shall contain 10% of their mass in the form of gas, 
and we assume that about 70% of the available gas in the disk is 
consumed in an exponentially declining starbursting phase. We use 
an e-folding time of 0.1 Gyr and an initi al star formation rate o f 
20 M0yr _1 to parametrize the starburst (Kazantzidis et al. 2005). 
We find that 300 Myr after the beginning of the starburst the 
colors of the m erger remnant are re d enough to fall into the selec- 
tion window of [van Dokkum (2005|). The elliptical galaxy is mod- 
elled with a 9 Gyr o ld single stellar population with red horizontal 
branch morphology (Maraston 2005). Both the old population and 
the starbursting population are assumed to have Salpeter initial 
mass functions and solar metallicity. All spectral energy d istribu- 
tions are derived from synthetic models (Maraston 2005). Mag- 
nitud es and colors are calculated with ZEBRA (Feldmann et al. 

2006) and normalized to Vega. Filter bands correspond to the 
Johnson-Cousins system. 

color difference between tidal debris and main elliptical's 
stars, and its standard deviation, are only of order ~ 0.3 
mag. The tidal tails produced by the accretion of disk 
galaxies are thus "red" in B — R color. 

Disk satellites responsible for the observed tidal fea- 
tures would likely host a gas component, which is 
not included in our ./V-body simulations. Such a gas 
component might cause a burst of star formation in 
the e arly phases of the interaction (jBarnes fc Hernquistl 
1991). Simulations show that tidally triggered stellar 
bars can develop soon after the first pericenter passage 
and drive substantial gas in flow within the central ~1 
kpc of the satellite ga laxy (jBarnes fc Hernquistl 119961 : 
Kazantzidis et al.l |2()()5i , where the gas is quickly ex- 
hausted in a burst of star formation. Star formation 
is also expected in the tidal tails at these early stages 
of the merger. Ram pressure stripping in the hot halo 
of the massive elliptical will however remove the reser- 
voir of diffu se gas and thus quench s tar formation at 
later times (jKawata fc Mulchaevi 120081 ) . Ram pressure 
becomes increasingly efficient as tidal shocks decrease the 
potential well of the satellite at each subsequent pericen- 
ter and is expected to remove most if no t all the gas 
that has not plunged into its central region (|Maver et al.l 

2007) . This is in agreement with the absence of a young 
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Fig. 10. — (Relative) number densities of disk and early-type galaxies, (a) Stellar mass functions of early-type (red dashed line) and disk 
(blue solid line) galaxies in the local universe ([Driver ct al. 2006). The vertical, black dot-dashed line corresponds to galaxies of 2.7 X 10 11 
Mq. The inset shows the ratio of the two stellar mass functions and indicates its value for a 4 times (magenta dashed line) and 10 times 
(green dot-dashed) less massive satellite galaxy. The ratio lies in the range of unity for the considered masses and increases with decreasing 
stellar mass, (ft) Relative number density fn{R) of disk vs. early- type galaxies for a given mass ratio R after integrating over all masses 
of the elliptical host galaxy (black solid line). Mass ratios of 1:4 and 1:10 are indicated by the green dot-dashed and the magenta dashed 
lines. The respective values of the average relative number density are 2.0 and 3.0, respectively. Under the assumption that the merger 
probabilities are proportional to the number densities, a typical 1:4 or 1:10 merger with an elliptical galaxy will involve much more likely 
a disk than another early-type system. 



stellar component in - and the red colors of - the ob- 
served tidal features. We predict, however, that very 
faint and diffuse gaseous tidal tails and plumes should 
be detectable around elliptical galaxies with red stellar 
tidal signatures. 

After the merger the newly formed stars in a disk- 
dominated progenitor might affect the overall B — R color 
of the remnant. Making a sensible guess for the param- 
eters of the central starburst, an e-folding time of 100 
Myr and an initial star formation rate of 20 M0yr _1 , we 
find that within 500 Myr the B — R color of the remnant 
will closely resemble the original B — R color of the el- 
liptical merger progenitor to less than 0.1 mag (Fig. [5]). 
In addition, a nonnegligible fraction of nearby elliptical 
galaxies are found to host a fe w, up to a few tens, p ercent 
of stellar mass in young stars (jTrager et al.l l2000). 

The existence of red tidal features around the most 



massive elliptical galaxies can thus be well explained with 
minor accretion events of disk-dominated satellite galax- 
ies. Consequently, our work supports an early (z > 1) 
mass assembly epoch for the most massive galaxies in 
the Universe with a later accretion of smaller mass units 
that do not contribute significantly to the total mass. 
Future observations will have to clarify the relative im- 
portance of dry and wet mergers at such early epochs in 
the formation of the bright elliptical galaxy population. 
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